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I  1  3-  ABSTRACT 


Measurements  of  low- temperature  thermal  conductivi ty  K  have  been  made  on  high- 
purity  Ge  single-crystal  samples  with  various  degress  of  surface  polish,  to  complement 
similar  measurements  previously  reported  [2]  on  silicon.  As  a  result  of  the  more  pro¬ 
nounced  isotopic  impurity  content,  bulx  scattering  could  not  be  completely  eliminated 
even  in  our  thinnest  (~0.5  mm)  and  roughest  samples.  Accordingly,  a  Calloway-type 
1-10)  analysis  was  used  in  order  to  extract  the  boundary-limited  mean  free  path. 

On  progressive  polishing,  the  conductivity  increased  markedly  and  deviated  markedly 
from  the  (approximate)  T3  behavior  of  the  rough  samples,  as  with  silicon.  A  good  fit 
to  the  measured  K  vs,  T  curves  was  obtained  by  taking  into  account  the  scattering  by 
the  isotopes,  the  polished  side  surfaces,  and  the  rough  end  surfaces.  For  the  side 
surfaces,  Ziman’s  theory  [12]  of  diffraction  by  the  surface  irregularities  in  an 
infinitely  long  crystal  was  employed.  Although  the  form  of  the  resulting  curves  fit 
the  data  well,  the  required  values  of  the  roughness  parameter  seem  anomalously  small. 
However,  the  silicon  data  are  equally  well  fitted  by  comparable  values  of  this  para¬ 
meter. 
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BOUNDARY  SCATTERING  OP  PHONONS  IN  GERMANIUM  AND  SILICON 

D.  R.  FrankI  and  George  J.  Campisi 
Department  of  Physics,  Pennsylvania  State  University 
University  Park,  Pennsylvania  16802  U.S.A. 

ABSTRACT 

Measurements  of  low-temperature  thermal  conductivity  ic  have  been  made  on 
high-purity  Ge  single-crystal  samples  with  various  degrees  of  surface  polish, 
to  complement  similar  measurements  previously  reported  [2]  on  silicon.  As  a 
result  of  the  more  pronounced  Isotopic  impurity  content,  bulk  scattering  could 
not  be  completely  eliminated  even  in  our  thinnest  (-0.5  mm)  and  roughest 
samples.  Accordingly,  a  Calloway-type  [10] analysis  was  used  in  order  to 
extract  the  boundary-limited  mean  free  path. 

On  progressive  polishing,  the  conductivity  increased  markedly  and 
deviated  markedly  from  the  (approximate)  T*  behavior  of  the  rough  samples,  as 
with  silicon.  A  good  fit  to  the  measured  ic  vs.  T  curves  was  obtained  by 
taking  into  account  the  scattering  by  the  isotopes,  the  polished  side 
surfaces,  and  the  rough  end  surfaces.  For  the  side  surfaces,  Ziman's  theory 
[12]  of  diffraction  by  the  surface  irregularities  in  an  infinitely  long 
crystal  was  employed.  Although  the  form  of  the  resulting  curves  fit  the  data 
well,  the  required  values  of  the  roughness  parameter  seem  anomalously  small. 
However,  the  silicon  data  are  equally  well  fitted  by  comparable  values  of  this 
parameter. 
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I.  INTRODUCTION 


At  low  temperatures  in  electrically  insulating  crystals,  only  two  classes 
of  phonon  scattering  processes  are  present:  "normal"  (momentum-conserving) 
three-phonon  processes,  and  collisions  with  fixed  imperfections.  Among  the 
latter  are  the  surface  of  the  crystal  and,  in  the  bulk,  dislocations,  point 
defects,  chemical  and  isotopic  impurities,  etc.  The  analysis  of  heat 
transport  in  the  presence  of  more  than  one  scattering  mechanism  is 
exceedingly  complicated.  Fortunately,  many  materials  are  now  obtainable  in 
sufficient  purity  and  perfection  that  the  bulk  defect  scattering  is 
negligible.  The  surface  is  then  the  primary  barrier  to  the  unimpeded  flow  of 
heat. 

In  recent  years,  there  has  been  a  marked  increase  in  interest  in  this 
"boundary-scattering"  regime,  with  regard  to  both  the  understanding  of  the 
fundamental  scattering  processes,  and  the  use  of  the  phenomena  as  a  probe  for 
the  surface  region  of  a  crystal.  Recent  experimental  measurements  include 
those  of  Thacher  [1]  on  LiF,  Hurst  and  FrankI  [2]  on  Si,  McCurdy  et  al  [3]  on 
CaF  and  Si,  and  Bird  and  Pearlman  [4]  on  Ge.  Most  of  these  measurements  were 
made  on  rough-surfaced  crystals.  With  the  exception  of  Reference  [2],  the 
early  work  of  Berman,  Foster,  and  Ziraan  [5],  and  the  present  work,  little 
attention  has  been  paid  to  the  extremely  interesting  case  of  highly  specular 
surfaces . 

The  theory  for  the  boundary-scattering  regime  is  rather  more  complicated 
than  ordinary  bulk  transport  theory,  since  not  only  the  kineraatical  properties 
of  the  carriers  but  also  the  points  of  origination  of  their  "flight  paths" 
must  be  taken  into  account.  Thus,  in  the  general  case,  one  is  faced  with  six- 
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dimensional  integrals  over  both  wave-vector  space  and  real  space.  However,  in 
the  absence  of  bulk  scattering,  all  flight  paths  originate  at  the  surface,  and 
the  integrals  reduce  to  five  dimensions.  These  have  been  evaluated  in  certain 
cases: 

1.  Casimir  [6]  considered  infinitely  long  circular  and  square  cylinders 
of  elastically  isotropic  material  with  perfectly  rough  (diffuse)  surfaces. 

His  classic  result  may  be  put  in  the  form 


Kiso 


■  |  c  <s>  *c. 


(1) 


where  C  is  the  specific  heat,  <S>  is  an  appropriate  average  sound  speed 
(averaged  over  the  three  acoustic  modes)  and  lc  is  the  so-called  "Casimir 
length",  given  by 

,  m  r2R  for  a  circular  cylinder,  ,2) 

c  i 1.12  W  for  a  square  cylinder,  v 

where  R  and  W  are,  respectively,  the  radius  of  the  circle  or  the  side  of  the 
square. 


2.  Berman,  Simon,  and  Ziman  [7]  extended  the  calculations  to  Include 
finite  length  L  and  non-zero  "specularity  factor"  P  for  circular  cylinders. 
The  result  was  that  fcc  was  replaced  by  a  corrected  length  l: 

l  -  &c  F(R/L,  P).  (3) 

For  an  infinitely  long  circular  cylinder,  the  correction  factor  is 

F(0,  P)  -  (1  +  P)/(l  -  P).  (4) 


f  3.  The  experimental  results  of  Hurst  and  Frankl  [2]  and  of  McCurdy  et  al 

[3]  as  well  as  certain  heat-pulse  effects  [83,  made  it  clear  that  the  elastic 
anisotropy  played  a  highly  significant  role.  The  key  factor,  as  pointed  out 
h  by  McCurdy  et  al,  is  that  the  phonons  travel  with  the  group  velocity  rather 

than  the  wave  velocity  (ordinary  sound  velocity).  On  this  basis,  they  were 
able  to  calculate  the  conductivity  of  finite-length,  rough-surfaced,  circular, 
square,  and  rectangular  cylinders  in  good  agreement  with  experiment.  Needless 
to  say,  their  results  reduced  to  Eq.  (1)  in  the  appropriate  limits. 


II.  EXPERIMENTAL 


A.  Method 


Oriented  germanium  rods  were  cut  from  a  single  crystal  of  high  chemical 
purity  and  of  the  natural  isotopic  mixture,  with  an  abrasive-wire  saw.  Rough 
surfaces  were  prepared  by  lapping  with  600-grit  SIC  (particle  size  ~9u)  on  a 
glass  block,  followed  by  sand-blasting.  Diamond  polishing  was  done  with  6|i 
and  lu  pt.wder  paste  starting  on  a  power-driven  cloth  lap  and  finishing  by 
hand  on  a  lens-paper  lap.  Further  polishing  was  done  with  a  commercial 
preparation  called  Lustrox.1'  Samples  were  square  cylinders,  typically  about 
1  mm  wide  x  40  mm  long,  and  usually  oriented  in  the  [100]  direction  with 
faces  of  {110}-type.  Thermal  conductivity  measurements  were  made  by  the 
standard  two-heater  method  in  a  He"  cryostat,  using  germanium  resistance 
thermometers.  Details  are  given  in  Reference  [2]. 


■^Sold  by  Tizon  Chemical  Corp., 


Flemington,  New  Jersey,  U.S.A. 
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B.  Results  for  Silicon 


The  results  previously  reported  [2]  may  be  summarized  as  follows: 

1.  For  rough-surfaced  crystals,  the  expected  ^-dependence  was 
reasonably  well  born  out;  actual  values  of  d(log  <)/d(log  T)  lay  mainly 
between  2.90  and  3.00.  The  size  dependence  followed  roughly  the  form  of 
Eq.  (2). 


2.  There  was  a  pronounced  anisotropy  in  the  conductivity,  the  [100] 
direction  being  roughly  twice  as  large  as  the  [110]  or  [111]. 

3.  On  polishing,  the  conductivity  increased  markedly,  especially  at  the 
lower  temperatures.  Apparent  P  values  >  0.9  were  obtained  in  a  few  eases. 


C.  Results  for  Germanium 


Figure  1  shows  typical  data  for  rough-surfaced  (sand-blasted)  rods.  It 
will  be  noted  that  d(log  ie)/d(log  T) falls  considerably  below  3.0  even  at  the 
lowest  temperatures.  This,  it  will  be  seen  later  is  a  result  of  the  isotope 
scattering.  The  size  effect  data  tend  to  substantiate  Eq.  (3)  within  about 
±5JS,  the  shape  correction  factors  being  close  to  1.0  for  these  samples.  The 
values  agree  well  (at  2°).  with  the  pure  [100]  sample  (sample  2)  of  Carruthers 
et  al  [9j.  In  contrast,  the  data  reported  by  Bird  and  Pearlman  [4]  for  their 
*pure"  sample  lie  about  20X  lower,  after  correcting  for  the  size  difference 
and  the  shape  factor  F  *  0.85.  Possibly  the  discrepancy  is  due  to  the  rather 
short  rods  used  by  these  authors. 

Figure  1  also  illustrates  the  anisotropy,  which  is  seen  to  be  somewhat 
less  pronounced  than  in  silicon.  In  passing,  it  may  be  mentioned  that  we  have 
also  observed  anisotropy  In  the  same  sense  in  LiF,  although  the  reproducibil¬ 
ity  of  the  data  was  not  very  good. 


Figure  1.  Thermal  conductivity  of 
rough-surfaced  germanium 
rods. 


Figure  2.  Effect  of  polishing  on 
thermal  conductivity  of 
germanium  [100]  rods. 
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The  effects  of  progressive  polishing  are  illustrated  in  Figure  2. 
Qualitatively  they  are  the  same  as  in  silicon,  although  the  more  serious 
isotope  scattering  lowers  the  conductivity  very  considerably.  It  will  be 
noted  that  the  enhancement  factor  is  again  considerably  greater  at  the  lower 
temperatures,  showing  that  the  specularity  factor  P  must  be  dependent  on  the 
phonon  wavelength. 


Several  attempts  '  %ve  been  made  to  improve  the  surface  smoothness  still 
further  by  chemical  etching  between  the  various  mechanical  abrasion  processes. 
The  intent  here  was  to  reduce  to  a  minimum  the  amount  of  damaged  material  in 
the  sample.  The  main  result  has  been  to  make  it  extremely  difficult  to  obtain 
reliable  data,  owing  to  very  poor  thermal  contact  of  the  thermometers  to  the 
samples.  The  same  difficulty  was  encountered  with  silicon,  but  not  as 
severely.  One  remedy  attempted  was  to  sandblast  a  small  strip  under  each 
thermometer  clamp.  However,  the  small  region  of  diffuse  surface  produced  an 
appreciable  decrease  in  conductivity,  so  this  approach  was  abandoned.  Work 
along  these  lines  is  still  in  progress. 


III.  ANALYSIS  OF  THE  DATA 


A.  Wavelength  Dependence  of  Specularity  Factor 


We  first  consider  an  approximate  analysis  in  which  the  anisotropy  effects 
are  Ignored.  The  main  features  dealt  with  are  (a)  the  increasing  departure 
from  Independence  as  the  surface  smoothness  increases,  and  (b)  the  separation 
of  surface  and  bulk  scattering.  The  standard  method  for  analysis  of  mixed 
scattering  is  via  the  theory  of  Callaway  CIO]  as  amended  by  Holland  [11]. 
Accordingly,  we  write 


k'T’/h1  p 


x*  ex  dx 


1*1  6ir2S 


o  (ex  -  1)2t“* 


(5) 


where  i  =*  1  designates  the  longitudinal  modes,  i  *  2,  3  the  transverse  modes, 
S.  is  an  appropriate  direction-averaged  sound  velocity  for  the  1th  mode-type, 
and  the  other  symbols  have  their  usual  meanings.  The  inverse  relaxation  times 
were  written  as 


<Ti>"‘  "  ^isotope^  +  +  <TB,i> 


-i 


(6) 


The  isotope  scattering  and  normal-process  scattering  terms  were  taken  from 
Holland  [11]  as 


^Tisotope^  *  Aw%  =  A(kBT/Ii)*x*, 


(6) 


and 


(THr*  *  Bu2T3  =  B(kB/h)2T5x2, 


,.2ml  _ 


2m5_2 


(7) 


with  A  *  1.32  x  10~  *  sec1,  B  =  3.8  x  10~2*  sec  deg~  for  silicon,  A  =  2.57  x 
10"**  sec*,  B  *  2.77  x  10"43  sec  deg"3  for  germanium.  The  values  of  (tn)"3 
were  negligible  throughout,  and  thus  the  exact  form  chosen  in  Eg.  (7) 
is  inconsequential. 


For  the  boundary  scattering,  we  wish  to  take  into  account  the  specularity 
of  the  surfaces  and  the  finite  length  of  the  rods.  As  noted  previously,  the 
progressive  departure  from  T3-dependence  makes  it  clear  that  the  specularity 
is  wavelength-dependent.  On  general  grounds,  it  could  be  expected  to  be 
angle-dependent  also,  but  there  is  no  way  of  incorporating  this  into  a  theory 
of  the  present  type.  Ziman  [12]  and  Soffer  [13]  have  considered  the  problem 
of  reflection  of  waves  from  rough  surfaces.  Soffer  obtains,  in  the  limit  of 
no  lateral  correlation  in  the  surface  roughness 


P(k)  *  exp  [-(2kneff)2], 


(8) 


% 


■fj 
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where  neff  =  n  cos  0lnc,  (8a) 

k  =  magnitude  of  phonon  wave-vector,  0.  *  angle  of  incidence,  and 

n  =  rms  height  deviation  in  the  surface.  Then,  in  view  of  Eq.  (3),  we 
take  as  a  plausible  averaged  expression 

=  si  [Wu  "  p)/(1  +  p)  +  (9) 

The  foregoing  expressions  were  gathered  into  Eq.  (5)  and  the  integrals 
computed  for  various  values  of  neff*  Results  are  shown  in  Figures  3  and 
In  these  calculations,  the  velocities  used  for  the  transverse  modes  were 
averages  of  the  type  calculated  by  Hasegawa  [14],  viz. ,  5.42  x  10s  cm/sec  for 
Si  and  3.29  x  10s  cm/sec  for  Ge.  There  is,  of  course,  no  reason  to  assume 
this  to  be  a  correct  average  in  the  present  context,  but  it  is  seen  that  the 
resulting  fit  for  the  rough-surfaced  samples  of  both  materials  is  reasonably 
good,  with  no  adjustable  parameters.  Clearly,  some  slight  adjustments  of  the 
velocities  could  improve  the  fit3. 


12  3  4  12  3  4 

T  (K)  TOO 


Figure  3.  Fit  of  Eqs.  (5)  -  (9)  to 

data  on  silicon.  Lowest  set 
of  points  for  rough  surface, 
next  for  ly  diamond  polish, 
upper  three  for  Lustrox 
polish. 


Figure  4.  Fit  of  Eqs.  (5)  -  (9)  to 
data  on  germanium.  Lowest 
set  of  points  for  rough 
surface,  next  for  ly 
diamond  polish,  upper  two 
for  Lustrox  polish. 


The  curves  for  partly  specular  surface  (finite  values  of  neff)  form  a 
family  that  seems  to  represent  the  observed  temperature  dependences 
reasonably  well.  Furthermore,  the  apparent  values  of  neff  are  comparable  for 
the  two  materials  at  each  stage  of  polishing-.  These  values  are,  to  be  sure, 
surprisingly  small  in  relation  to  the  grit  sizes  used  in  the  polishing.  They 
can  possibly  be  understood,  however,  as  a  result  of  the  cos  0ipC  factor  in 
Eq.  (8a),  which  would  lead  to  an  overwhelming  preference  for  those  phonons 
traveling  nearly  parallel  to  the  rod  axis. 


Directional  Dependence  of  Specularity  Factor 


In  order  to  help  decide  whether  the  basic  mechanism  (i.e.,  diffraction) 
discussed  in  the  preceding  section  may  be  the  correct  one,  it  would  be 


.  •#* 


*** 
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desirable  to  extend  the  anisotropic  transport  theory  of  McCurdy  e£  al  [3]  to 
the  case  of  partly  specular  reflection.  However,  the  problem  is  a  very 
complicated  one.  Among  the  effects  to  be  considered  are  the  following: 

1.  The  factor  (1  +  P)/(l  -  P)  is  valid  only  for  an  isotropic  circular 
rod,  in  which  successive  reflections  yield  identical  flight  paths.  In  a  square 
rod,  there  are  different  reflected  flight  paths  depending  on  the  original 
point  of  emission,  and  there  are  two  distinct  angles  of  incidence,  hence  two 
P-values . 

2.  The  result  of  a  specular  reflection  is  to  preserve  momentum 
components,  not  velocity  compoenents.  A  consequence,  according  to  Price  [15], 
is  that  specular  reflection  does  produce  resistance.  Thus,  the  (1  +  P)/(l  -  P) 
factor  is  wrong  even  in  the  limit  P  -*■  1. 

3.  When  a  single  mode  is  reflected,  the  general  result  is  that  at  least 
two  reflected  waves  are  required  in  order  to  satisfy  the  boundary  condition 
[16] .  Thus,  instead  of  a  single  function  P(k),  we  need  a  matrix  Pjj(k) 
connecting  the  various  modes. 

‘’.‘onceivably,  this  last  point,  could  hold  the  key  to  the  explanation  of 
our  anomalous  values  of  neff-  In  effect,  we  find  that  the  polishing  increases 
the  conductivity  by  an  unexpectedly  large  factor.  This  could  occur  if  there 
were  a  net  transfer  of  phonons  to  more  "efficient”  modes.  In  fact,  according 
to  Thurston  (16],  this  is  Just  what  occurs:  a  transverse  wave  (slow), 
produces  a  reflected  longitudinal  wave  (fast),  which  moves  more  nearly 
parallel  to  the  surface.  Beyond  a  critical  angle,  the  latter  becomes  a 
Rayleigh  wave,  which  might  in  a  good  sample  be  virtually  un-attenuated . 
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